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ABSTRACT 

Strong lensing has developed into an important astrophysical tool for probing 
both cosmology and galaxies (their structure, formation, and evolution). Using 
the gravitational lensing theory and cluster mass distribution model, we try to 
collect a relatively complete observational data concerning the Hubble constant 
independent ratio between two angular diameter distances Dd s /D s from various 
large systematic gravitational lens surveys and lensing by galaxy clusters com- 
bined with X-ray observations, and check the possibility to use it in the future as 
complementary to other cosmological probes. On one hand, strongly gravitation- 
ally lensed quasar-galaxy systems create such a new opportunity by combining 
stellar kinematics (central velocity dispersion measurements) with lensing geom- 
etry (Einstein radius determination from position of images). We apply such a 
method to a combined gravitational lens data set including 70 data points from 
Sloan Lens ACS (SLACS) and Lens Structure and Dynamics survey (LSD). On 
the other hand, a new sample of 10 lensing galaxy clusters with redshifts ranging 
from 0.1 to 0.6 carefully selected from strong gravitational lensing systems with 
both X-ray satellite observations and optical giant luminous arcs, is also used 
to constrain three dark energy models (ACDM, constant w and CPL) under a 
flat universe assumption. For the full sample (n = 80) and the restricted sample 
(n = 46) including 36 two-image lenses and 10 strong lensing arcs, we obtain 
relatively good fitting values of basic cosmological parameters, which generally 
agree with the results already known in the literature. This results encourages 
further development of this method and its use on larger samples obtained in the 
future. 
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Subject headings: Gravitational lensing: strong - (Cosmology:) cosmological pa- 
rameters - (Cosmology:) dark energy 



Introduction 



Pioneering observations of type la supernovae (SNe la) (IRiess et al. Il998t iPerlmutter et al.l 



19991 ) have demonstrated that our present universe is passing through an accelerated phase 
of expansion preceded by a period of deceleration. A new type of matter with negative 
pressure known as dark energy, has come up to explain the present phase of acceleration. 
The simplest candidate of dark energy, the cosmological constant (A), is consistent with var- 



ious observations suc h as more precise supern ova data (jRiess et al. 



Kowalski et al.ll2008l ). the CMB observations (ISpergel et al. 



2007 



light elements abundance from Big Bang Nucleosynthe sis (IBurles et al. 



aco ustic oscillations ( BAO) detected in SDSS sk y survey (lEisenstein et al. 



ies (iDaly et al.l 120091 ) . and gamma-ray bursts (lAmati et al 



2004 



Davis et al. 



Komatsu et al. 



2007 



2009J), the 



2001), the baryon 



200 



20051 ) , radio galax- 



However, various other 
models were proposed a s candidates of dark energy, such as the t ypical dynamic al scalar 
field called quintes sence (|Caldwell et al.lll998l) . phanto m corrections (Caldwell 2002 ), a joint 



quintom scenario (IFeng et al 



Biesiada. Godlowski fe Szydlowski 



2005) or Chaplygin gas ( Kamenshchik et al 



2001 



Zhu 2004 



2005t IZhang feZhul 120061 ). to mention just a few out of a 



long list. On the other hand there are still many ot her ways to understand the accelerat- 
ing universe, such as Modified Friedmann Equati on (IFreese &; Lewis! 120021 : IZhu et al.l 120041 ) 
and Dvali-Gabadadze-Porrati(DGP) mechanism (jDvali et al.l 120001 ) . But until now none of 
these models was demonstrated superior over the other. Besides, while updating the current 
estimates of cosmological model parameters, one should try to use new probes. Strongly 
gravitationally lensed systems belong to this category. They can provide the information 
on two angular diameter distances, Dd s and D s . One is the distance to the source and the 
other is that between the defector and the source. Since angular diameter distance depends 
on cosmological geometry, we can use their ratios to constrain cosmological models. 



The discovery of strong gravitational lensing in Q0957+561 (I Walsh et al.l 1 19791 ) opened 
up an interesting possibility to use strong lens systems in the study of cosmology and as- 
trophysics. Up to now, stro ng lensing ha s developed into an important astrophysical tool 



for probing both cosmology (IZhu 



Zhu fc Mauro 



2001 



2008a 



tion) (IZhu fc Wu 



Ofek et al 



Zhu et al. 



1997 



2003 



2000a b; Chae 2003; Chae et al. 



20041 : Mitchell et al.ll2005 



2008bh and galaxies (their structure, formation, and evolu- 



Mao fc Schneider 1998; Jin et al. 2000; Keeton 2001; Kochanek fc White 



Treu et al.ll2006a[ ). Now several hundreds of strong lens systems pro- 
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duced by massive galaxies have been discovered, but only ~ 90 galactic-scale strong lenses 
with known redshift of the lens and the source and measured image separation can form 
well-defined samples useful for statistical analyses. These well-defined strong lenses are 
particularly useful not only for constraining the statistical properties of galaxies su ch as 
stellar velocity dispersions or galaxy evolution (IChae fc Mad 120031 ; lOfek et al.l 120031 ). but 
also for constraining cosmological parameters such a s the present-day matter densit y Q m , 
dark energy density Q x and its equation of state w ( IChad 120031 ; Mitchell et al.ll2005l ). For 
example, the Cosmic Lens All-Sky Survey (CLASS) statisti cal data, which consists of 8958 
radio sources out of which 13 sources are multiply imaged (IBrowne et al.ll2003l ; IChad 120031 ) 
was first extensively used by IChad (120031). who found Q m ~ 0.3 assuming a flat cosmology 



and non-evolving galaxy populations. iMitchell et al.l (120051 ) reused this CLASS statistical 



sample based on the velocity dis persion function (VDF) o f early-type ga l axies d erived from 
the SDSS Data Release 1 (DR1; Istoughton et al.l (|2002h ). Izhu fc Maurol (l2008at ) reanalyzed 
10 CLASS multiply-imaged sources whose image-splittings are known to be caused by single 
early-type galaxies to check the validity of the DGP model with radio-selected gravitational 
lensing statistics. More recently, the distribution of gravitationally-lensed image separations 
observed in the Cosmic Lens All-Sky Survey (CLASS), the PMN-NVSS Extragalactic Lens 
Survey (PANELS) , the Sloan Digital Sky Survey (SDSS) and other surveys was used by 
Cao fc Zhul ( I2011af ). who found w < —0.52 assuming a flat cosmology and adopting semi- 
analytical modeling of galaxy formation. The idea of using strongly gravitationally lensed 
systems, in particular measurements of their Einstein radii combined with spectroscopic 
data, for measuring cosm ological parameters including the cosmic equation of state was dis- 
cussed in Biesiada (12006[) and also in a mo re recent papers (jGrillo et aL||2008j ; iBiesiada et al.l 
2010l : iBiesiada. Malec k Piorkowskalboilh . 



On the other hand, galaxy clusters, as the largest dynamical structures in the universe 
are also widely used both in cosmology and astrophysics. Firstly, their mass distributions 



at di fferent redshifts can be described by the Press-Schechter function (IPress &: Schechter 



1974 ). which reflects the linear growth rate of density perturbations and therefore can pro- 



yide constraints on cosmological parameters such as the matter and da rk energy densities 



(Bo rgani et al.lll999l ). Secondly, combining the Sunyaev-Zel'dovich effect (jSunyaev fc Zeldovich 



Hubble constant and other cosmological parameters in given cosmological model ( 


Reese et al. 


2002; 


Schmidt et al. 


2004; 


Jones et al. 


2005; 


Bonamente et al. 


2006; 


Zhu & Fujimoto 


2004 





Relevant discus sions on the corr e ctions to the Sunyaev- Zeld ovich effect for galaxy clusters 



can be found in lltoh et al.l (119981 ); iNozawa et al.l ( 119981 . 120061 ). More importantly, giant arcs 



generated by the galaxy cluster are perfect indicators of its surface mass density, while the 
mass distribution of the cluster's mass halo can be modelled from X-ray luminosity and 
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temperat ure, which may p rovide certain observable (jSerenol 120021 ; ISereno fc Longol I2004J ). 
Recently, lYu fc Zhul ( 120101 ) collected a new sample with such data from an online database 
BAX and various literatur e , whi ch led to some interesting results compared with those ob- 
tained by lSereno &: Longol (120041 ) . 



In this paper, we try to collect a relatively complete observational data concerning the 
Hubble constant independent ratio between two angular diameter distances D^ a jD s from 
various large systematic gravitational lens surveys and galaxy cluster data. This paper is 
organized as follows. In Section [2J we briefly describe the methodology for both strong grav- 
itationally lensed systems: galactic lenses and galaxy clusters. Then, in Section [3] we present 
the Dds/D s data from various large systematic gravitational lens surveys and lensing galaxy 
clusters with X-ray observations and optical giant luminous arcs. We further introduce three 
popular cosmological models tested in Section HI Finally, we show the results of constraining 
cosmological parameters using MCMC method and conclude in Section [5j 



2. The Method 

Gravitational lensing is one of the successful predictions of General Relativity. Strong 
gravitational lensing occurs whenever the source, the lens and the observer are so well 
aligned that the observer-s ource direction lies inside the so-called Einstein ring of the lens. 



Paczynski fc Gorskil ( 119811 ) tried to use lensing images as indicators to estimate cluster mass 



and constrain cosmological constant. 

In a cosmological context the source is usually a quasar with a galaxy acting as the lens. 
Strong lensing reveals itself as multiple images of the source, and the image separations in 
the system depend on angular diameter distances to the lens and to the source, which in 
turn are determined by background cosmology. Since the discovery of the first gravitational 
lens the number of strongly lensed systems increased to a hundred (in the CASTLES data 
base) a nd is steadily incre asing following new surveys like the Sloan Lens ACS (SLACS) 



survey (INewton et al.ll201ll ). This opens a possibility to constraining the cosmological model 



provided that we have good knowledge of the lens model. 

Now, the idea is that the formula for the Einstein radius in a SIS lens (or its SIE 
equivalent), 

D A (z,z s ) g 2 SIS 

depends on the cosmological model through the ratio of (angular diameter) distances between 
lens and source and between observer and lens. Under flat Friedman- Walker metric, the 
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angular diameter distance reads 



dz' 



H (l + z)J E(z';p) 



(2) 



where Hq is the Hubble constant and E(z; p) is a dimensionless expansion rate dependent 
on redshift z and cosmological model parameters p. If the Einstein radius 9e from image 
astrometry and stellar velocity dispersion o (or central velocity dispersion er ) from spec- 
troscopy can be determined, this method can be used to constrain cosmological parameters. 
The advantage of this method is that it is independent of the Hubble constant value and is 
not affected by dust absorption or source evolutions. However, it depends on the measure- 
ments of o"o and lens modelling (e.g. singular isothermal sphere (SIS) or singular isothermal 
ellipsoid (SIE) assumption). Hopefully, spectroscopic data for central parts of lens galaxies 
became available from the Lens Structure and Dynamics (LSD) survey and the more re- 



(Treu et al. 


2006a 


b; 


Grillo et al. 


2008) 



assumption in gravitational lensing studies and should be accurate enough as first-order ap- 
proximation to t h e mea n properties of galaxies relevant to statistical lensing. For example, 



Koopmans et al.l (120091 ) found that inside one effective radius massive elliptical galaxies are 



kinematically indistinguishable from an isothermal ellipsoid. In the previous wor ks, such an 



1996; 


King et al.l 


1997; 


Fassnacht & Cohen 


1998 


2006. 


2009; 


Treu et al. 


2006a 


b 


)• 



However, let us note here that the velocity dispersion a sis of the mass distribution 
and the observed stellar velocity dispersion o"o need not be the same. IWhite &; David (119961 ) 
argued that there is a strong indication that dark matter halos are dynamically hotter than 
the luminous stars based on X-ray observations, and dark matter must necessarily have a 
greater velocity dispersion than the visible stars. In this paper, we adopt a free parameter 
£e that relates the vel ocity dispersion asis and the stellar velocity dispersion <To (IKochanek 
19921 : lOfek et~al~1l2003h : 

vsis = fE&o- (3) 

To be more specific, we have kept f& as a free parameter, since it mimics the effects of: 
(i) systematic errors in the rms difference between o"o (observed stellar velocity dispersion) 
and asis (SIS model velocity dispersion); (ii) the rms error caused by assuming the SIS 
model in order to translate the observed image separation into 9e, since the observed image 
separation does not directly correspond to 9e] (hi) so ftened isothermal sphere pote ntials 
which tend to decrease the typical image separations (INarayan fc Bartelmannl Il996l ). and 
could be represented by Je somewhat smaller than 1. 
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Martel et al.l (120021 ) found that the presence of the background matter tends to increase 
the image separations produced by lensing galaxies from ra y-tracing s i mulat ions in CDM 
models, though this effect is small, of order of 20% or less. IChristleinl (120000 showed that 
richer environments of earl y type galaxie s may have a higher ratio of dwarf to giant galaxies 



than the field. However, iKeeton et al.l ( 120001 ) showed that this effect nearly cancels the 



effect of the background matter, making the distribution of image separations significantly 
independent of environment. They predicted that lenses in groups have a mean image 
separation which is ~ 0.2" smaller than that of lenses in the field. Therefore, all these above 
factors can possibly affect the images s eparation by up to ±20%, which may be mimicked 
by introducing (0.8) 1 / 2 < f E < (1.2) 1 / 2 (bfek et al.lhool l 



In the method we use, the cosmological model enters not through a distance measure 
directly, but rather through a distance ratio 



T> th (z d ,z s ;p) 



D ds _i:;[dz'/E( Z '- P )} 



D s f z °[dz>/E(z';p)) 
and respective observable counterpart reads 



(4) 



V 



oh* 



c 2 e f 



(5) 



with its corresponding uncertainty calculated through propagation equation concerning the 
errors both on the s tellar velocity disp ersion a and the Einstein radius 6 E (~ 5% error for 
the Einstein radius (jGrillo et al.ll2008f )). 



Another source of systematic errors in our method comes from the fact that Einstein 
radius estimation from observed image positions depends on the lens model (SIS or SIE or 
the other realistic mass distribution). Moreover, the image separation could be affected by 
nearby masses (satellites, neighbor galaxies) or the structures along the line of sight. This 
last issue will also be discussed in the last section. Here, let us note that formally, at the level 
of Eq.flSD the f E factor does the double job accounting for systematics associated both with 
do and 9 E . Since the main goal of this paper is to constrain cosmological parameters, we 
firstly consider f E as a free parameter, obtain its best-fit value and probability distribution 
function P(f E ), and then treat it as a "nuisance" parameter to determine constraints on the 
relevant cosmolo gical parame t ers of i nterest. The proced u re of marginalization is carried out 
following that of Allen et al. ( 2008 ): Samushia Sz Ratral ( 2008 ). where -P(/e) is normalized 
to one and is usually taken to be a Gaussian or a 6(f E ) function peaked at the best-fit value 
of f E . We then integrate the likelihood function, 



£(p) 



C(p,f E )P(f E )df E 



(6) 
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and determine the best-fit values and confidence level contours from £(p). 



Moreover, strong lensing by clusters with galaxies acting as sources can produces giant 
arcs around galaxy cluste rs, which can also be used to constrain clusters' projected mass and 
cosm ological parameters (ILynds fc Petrosianlll986l ; iBreimer fc Sanderdll992t ISereno fc Longo 
20041 ). When a galaxy cluster is relaxed enough, t he hydrostatic isothermal spherical sym- 
metric fi- model (ICavaliere fc Fusco-Femianolll976l ) can be used to describe the intracluster 
medium(ICM) density profile: n e (r) = n e0 (1 + r 2 /r 2 ) 3/?X//2 ; where n e0 is the central electron 
density, fix an d r c denote the slope and the core radius, respectively. Assuming that whole 
gas volume is isothermal (with the temperature Tx), the gravity of relaxed cluster and its 
gas pressure should balance each other according to the hydrostatic equilibrium condition. 
With the approximation of spherical symmetry, the cluster mass profile can be given by 
M(r) = 3k Q^^f x r 2+ T i , where k B , m p and [i = 0.6 are, respectively, th e Boltzmann constant, 
the proton mass, and the mean molecular weight (IRosati et al.l 120021 ) . A theoretical surface 
density can be derived as 



3 k B T x (3 x 1 



(7) 



Combining this with the critical surface mass density for lensing arcs ( [Schneider et al.lll992[ ) 



-'obs 



4nGD d D ds 



+ 1, 



a Hubble constant independent ratio can be obtained 



V 



obs 



D 



ds 



obs 



6n ksTx(3 



X 



(9) 



The X-ray data fitting results may provide us the above mentioned relevant parameters 
such as T x , fix, and 9 C . The position of tangential critical curve 6 t is usually deemed to 
be equal to the observational arc position 9 arc . In this paper we assume that the deflecting 
angle has a slight dif ference with the arc radius angle, 8 t = e8 arc , with the correction factor 
e = (l/vT2) ± 0.04 ( jOno et al.lll999l ). The complete set of standard pri ors and allowance s 
of the above parameters included in Eq. [9] can be found in Table 1 of lYu fc Zhul ( 120101 ) . 
The observational T> obs and its corresponding uncertainty are also calculated through Eq. 

We stress here that the observational distance ratio T> has both advantage and disad- 
vantage. The positive side is that the Hubble constant H gets cancelled, hence it does not 
introduce any uncertai nty to the results. T he disadvantage is that the power of estimating 



fl m is relatively poor (IBiesiada et al.l 120101 ) . Therefore we only attempt to fit Q m in the 
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case of a ACDM model (where it is the only free parameter in flat cosmology). Then for 
both cases above (Eq. j5] and Eq. j9]), we can fit theoretical models to observational data 
by minimizing the y 2 function 

2^_V-Pf(P)- I, f SV2 



X IP 



^ ' . (io) 



where the sum is over the sample and a^, i denotes the variance of T> 



3. Sample used 



For the Einstein ring data, we first use a combined sample of 70 strong lensing systems 
with go od spectroscopic measurements of centr al dispersions from the SLACS and LSD 



surveys (IBiesiada et al. 



2010 



be found in ( jKoopmans fc Treu 



2002 



2003 



Newton et al.ll201lh. Original data concerning the sample can 



Treu fc Koopmand 12004 ; iTreu et al.l l2006al ) (see 



for details). In our sample of 70 lenses some have 2 images and some have 4. There are some 
general ar guments in favor of SIS model, but strictly speaking SIS lens should have only 
2 images ( IBiesiada et al.l l2010t IBiesiada. Malec fc Pi6rkowskall201ll ). so one can try to use 
only 2 image systems out of the full sample. Therefore we selected a subsample of n = 36 
lenses, which is summarized in Table [1] where the names of lenses in the restricted sample 
are given in bold. 

As for the strong lensing arcs, redshifts and temperatures of the galaxy clusters are 
always searched out directly from online databases, such as CDS (Th e Strasbourg astr onom- 
ical Data Center) or NED (NASA/IPAC Extragalactic Database). lYu k Zhul f)201Clh have 
chosen a new database established especially for X-ray galaxy clusters - BAX, which pro- 
vides detailed information including and 8 C . They also used the fitting results of Chandra, 



ROSAT, ASCA satellites and VIMOS-IFU su rvey flOta fc Mitsudal 12004 ; Bonamente et al. 



20061 ; ICovone et al.ll2005l ; iRichard et al.l 120071 ). The final statistical sample satisfy the fol- 



lowing well-defined selection criteria. Firstly, the distance between the lens and the source 
should be always smaller than that between the arc source and the observer, D ds /D s < 1, 
which rules out half of selected lensing arcs. Secondly, the arcs whose positions a re too far 
from characteristic radius {6 arc > 39 c ) should also be discarded (jYu fc Zhull2010i ). At last 



Yu fc Zhul (120101 ) obtained a sample of 10 giant arcs with all necessary parameters listed in 
Table [U 

Now the observational Dd s /D s data containing 80 data points for cosmological fitting 
are summarized in Table [H with errors calculated with error propagation equation. We also 
list a restricted sample containing 46 data points, which consists 36 two-image lenses and 10 
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strong lensing arcs. 



4. Cosmological models tested 



All cosmological models we will consider in this paper are currently viable candidates 
to explain the observed acceleration. Given the current status of cosmological observations, 
there is no strong reason to go beyond the simple, standard cosmological model with zero 
curvature and cosmological constant A (except for the conceptual problems arising when 
one attemp ts to reconcile i ts observed value with some estimate derived from fundamental 
arguments (lWeinbergill989f )). However, it is still interesting to investigate alternative models. 
And we hope that future observations of more accurate Dd s /D s data could allow to better 
discriminate various competing candidates. In the MCMC simulations we assume for each 
class of models the best fit values of parameters found in the present work, and vary them 
within their 2a uncertainties. We assume spatial flatness of the Universe throughout the 
paper, since it is strongly supported by independent and precise experiments e.g. a combined 
5-yr Wilkinson Microwave Anisotropy Probe (WMAP5 ), baryon acoustic osc illations (BAO) 
and supernova data analysis gives VL to t = 1.0050lo!oo6i fainshaw et al.ll2009l ). Moreover, the 



Vt m = 0.27 prior is used except in the ACDM model where the fit is attempted. 

For comparison we also performed fits to the newly released Union2 SNe la data (n=557 
supernovae) from the Su pernova Cosmology project covering a redshift range 0.015 < z < 1.4 
(lAmanullah et al.ll2010l ). In the calculation of the likelih ood from SNe la, we marginalize 
over the nuisance parameter (iDi Pietro fc Claeskensll2003l ) 



XsNe 



a B 2 1 ( G 



:n) 



where A = £f 7 (> data - ^f/af , B = £f ' (// ata - /i th )/a,f, C = £f l/crf, and the 
distance modulus is \i — 5 log(<ik/Mpc)+25, with the la uncertainty <jj from the observations 
of SNe la; and the luminosity distance di as a function of redshift z 



di 



:i+z) 



cdz' 



o H E{z';p) • 



;i2) 



4.1. The standard cosmological model (ACDM) 



We start our analysis by first setting out the predictions for the current standard cosmo- 
logical model. In the simplest scenario, the dark energy is simply a cosmological constant, 
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A, i.e. a component with constant equation of state w = p/ p = —1. If flatness of the FRW 
metric is assumed, the Hubble parameter according to the Friedmann equation is 

E\z;p) = n m (l + zf + n A , (13) 

where Q m and fl\ parameterize the density of matter and cosmological constant, respectively. 
Moreover, in the zero-curvature case (Q = Q m +Q\ = 1), this model has only one independent 
parameter: p = fl m . 



4.2. Dark energy with constant equation of state (wCDM) 

Allowing for a deviation from the simple w — — 1 case, the accelerated expansion is 
obtained when w < —1/3. In a zero-curvature universe, the Hubble parameter for this 
generic dark energy component with density Q x then becomes 

E\z- p) = fi m (l + zf + 0,(1 + zf^ w \ (14) 

Obviously, when flatness and Q m = 0.27 are assumed, it is a one-parameter model with the 
model parameter: p = {w}. 



4.3. Dark energy with variable equation of state (CPL) 

If the equation of state of dark energy is allowed to vary with time, one has to choose 
a suitable functional form for w(z), which in g eneral involves certain parametrization. N ow, 
we consider the commonly used CPL model ( Chevalier fc Polarski 2001 : Linderi I2OO3I ) . in 
which the equation of state of dark energy is parameterized as w(z) = u>o + w a j^, where 
w and w a are constants. The corresponding E(z) can be expressed as 



e 2 {z- P ) = a m (i + zf + (1 - n m )(i + z 



.3(1+WQ+W a ) 



cxp 



3WgZ 

l + z 



(15) 



There are two independent model parameters in this model: p = {wq, w a }. 



5. Results and conclusions 

In the first case, we consider as a free parameter and show the constraint results with 
the full n = 70 and the restricted n = 36 two-image galaxy lenses in Fig.[T]and Fig.|2j In order 
to derive the probability distribution function for the cosmological parameters of interest, 
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Fig. 1. — The 68.3 and 95.8 % confidence regions for ACDM model in the (fi m ,/e) plane 
obtained from the full n = 70 and the restricted n = 36 two-image galaxy lenses. The crosses 
represent the best-fit points. 




Fig. 2. — The 68.3 and 95.8 % confidence regions for wCDM model in the (w,/e) plane 
obtained from the full n = 70 and the restricted n = 36 two-image galaxy lenses. The 
crosses represent the best-fit points. 
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Fig. 3. — The marginalized constraint on Q m of ACDM model from 80 full Dd s /D s data and 
46 restricted Dd s /D s data. 




Fig. 4. — The marginalized constraint on w of u>CDM model from 80 full Dd s /D s data and 
46 restricted D ds /D s data. 
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we marginalize through Eq. [6] and perform fits of different cosmological scenarios on 
the full n = 80 sample as well as the restricted n = 46 sample with the results displayed in 
Table H 

For the full n = 80 sample (containing 70 galaxy lenses and 10 strong lensing arcs), 
first, in ACDM model where Q m is the only free parameter we were able to make a reliable 



fit on the samples considered. This result is a considerable improvement over iBiesiada et al. 
( 120101 ). where the authors failed to constrain Q m with their sample of twenty Einstein rings. 



Let us compare our results with previously known ones. T he current best fit value from 
cosmological ob s ervat ions is: Q\ = 0.73 ±0.04 in the flat case ( jPavis et al.ll2007l ). Moreover, 
Komatsu et al.l (120091 ) gave the best-fit parameter Q m = 0.274 for the flat ACDM model 
from the WMAP5 results with the BAO and SNe la Union data. We find that our value of 
Q m (see Table 2) obtained from the D ds jD s data is consistent with the previous works at 
la. Secondly, the best fit for the wCDM parameter agrees with that inferred from SNe la or 
WMAP5, and the ACDM model (w = -1) still falls within the la interval from the D ds /D s 
sample. Hence the agreement is quite good. Thirdly, concerning the evolving equation of 
state in the CPL parametrization, confidence regions in the (wo,w a ) plane are shown in 
Fig. P One ca n see t hat fits for w and w a are greatly improved as compared with those of 
Biesiada et al. (hoioh. The values inferred are also in agreement with the WMAP5 results 
presented in iHinshaw et al.l ( 120091 ) including combined WMAP5, BAO and SNe la analysis. 
Moreover, it can be seen that the concordance model (ACDM) is still included at la level 
for the Dds/D s data applied her e. For comparison we a lso plot the likelihood contours with 
the Union2 SNe la compilation (lAmanullah et al.ll2010l ). One can see that the w coefficients 
obtained from the Dd s /D s sample agrees with the respective values derived from supernovae 
data (almost the whole 2a confidence interval for w from the Union2 data set lies within 
the 2a CI from the D^/ 'D s data), which demonstrates the compatibility b etwee n the SNe la 
and Dds/D s data. This is also a great improvement over IBiesiada et al.l ( 120101 ). where SNe 
la results and strong lensing results were found marginally inconsistent at 2a. 

Working on the restricted n = 46 sample (containing 36 two-image lenses and 10 strong 
lensing arcs), despite the sample size has decreased dramatically, we find that fits on Q m in 
ACDM model are consistent with the standard knowledge (see Fig. [3]) and the best fit for 
the w parameter in quintessence scenario is higher than inferred from SNe la or WMAP5 
(see Fig. H]). Moreover, for the fits on wo and w a in CPL parametrization, even though 
confidence regions get larger in Fig. [5j the result also turns out to agree with SNe la fits. 
One should also note, that a systematic shift downwards in the (w ,w a ) plane persists. Such 
a shift in best-fitting parameters inferred from supernovae (standard candles, sensitive to 
luminosity distance) and BAO (stand ard rulers, sensitive to an g ular diameter distanc e) has 
already been noticed and discussed in lLinder fc Roberts! (120081 ); IBiesiada et al.l (120101 ). Our 
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result suggests the need for taking a closer look at the compatibility of results derived by 
using angular diameter distances and luminosity distances, respectively Recent discussions 
on the ideas of testing the Et h eringt o n reciprocity relati o n between these two distances ca n 
be found in Bassett fc Kunzl (120041): lUzan et all JpOoi h Eolanda. Lima fc Ribeirol feoid ): 
Cao fc Z hu J2011bh : [piorkowska. Biesiada fc Maled fcoioh . 



In conclusion our resul t s demonstrate that the method extensively investigated in lBiesiada 



(|2QQa); iGrillo et all (120081 ): Biesiada et all (l2010h : lYu fc Zhul (12010( 1 on simulated and ob- 
servational data can practically be used to constrain cosmological models. Moreover, good 
quality measurements of the relevant observational qualities such as the velocity disper- 
sion and Einstein radius turn out to be crucial. Finally, four important effects, neglected 
here, should be mentioned. One is that both the Einstein rings and X-ray observations of 
our new lensing sample come from different surveys or satellites (SLACS, LSD and SBAS 
and Chandra, ROSAT and ASCA, respectively), the differences in detectors and observing 
strategies may cause systematical errors which are hard to estimate. The second is that 
the observed image separation is affected by secondary lenses (satellites, nearby galaxies, 
groups, etc) in many cases. In this case, those lenses should not be used or the true $e 
corresponding to <Tq should be estimated through realistic modelling. However, most of our 
samples come fro m the SLACS survey where the role of environment has been assessed in 
Treu et all (120091 ). Namely, it was found that for SLACS lenses the typical contribution 
from external mass distribution is no more than a few percent. The third important effect 
is, that the statistical procedure for cluster lenses relies on many simplifying assumptions. 
The realistic errors should be estimated by more realistic model of galaxy clusters besides 
the hydrostatic isothermal spherical symmetric /3-model.The last one is the influence of line- 
of-sight mass contamination, with the signif icant effect of the large-scale structure on strong 

997). More recent results on this issue can be found 



in 



le nsing (IBar-Kana 



1996 



Keeton et al 



Dalai et all ( bood ): iMomcheva et al 



(120061 ) . In this paper, large scale structure effects 



which change the typical separation between images are also included in the parameter /g 

i ii 

an increase of an arbitrary order f% in the ve locity dispersion is equivalent to an increase 



of fs in the typical separation 9 (i.e., 9 oc a 2 ) (jMartel et al. 2002 ). In order to be complete 
with the discussion of possible errors one should also notice that the redshifts z s and zi 
are also known with some accuracy 5z s and 5zi which propagates into theoretical distance 
ratio calculations. In principle one should have accounted for them by suitable num erical 
simulations. However, based on the experience gained on SNIa ( iPerlmutter et al.lll999l ). this 
effect is likely to be much smaller than systematic errors discussed above. Another straight- 
forward solution based on Poissonian statistics suggests that a sample size of order of a few 
hundred lenses m ight reduce the line-of-sight 'noise' contamination down to a few percent 
( IKubo et al.ll2010l ). However, our Da s /D s data set is really small, and its range of redshift is 
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also limited. Fortunately, with the ongoing of various systematic gravitational surveys and 
more giant arc sur vey projects carried out by the International X-ray Observatory (IXO) 
( White et al.ll2010[). extended Roentgen Survey with an Imaging T elescope Array (eRosita) 
jPredehl et al.lboioh and the Wide Field X-ray Telescope (WFXT) flMurrav fc WFXT Team 
201oh being under way, the sample of strong lenses is growing rapidly, which may ease the 



problem of line-of-sight contamination. Future observations will definitely enlarge our set 
and make the method applied in this paper more powerful. 
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Table 1: Values of T> = D c i s /D s from lensing galaxy clusters and combined SLACS+LSD 
lens samples. The two-image lenses are written in bold. 



Systcm 



V 



rcf 



MS 0451.6-0305 0.550 2.91 0.785 0.087 

3C220.1 0.61 1.49 0.611 0.530 

CL0024.0 0.391 1.675 0.919 0.430 

Abell 2390 0.228 4.05 0.737 0.053 

Aboil 2667 0.226 1.034 0.837 0.124 

Abell 68 0.255 1.6 0.982 0.225 

MS 1512.4 0.372 2.72 0.734 0.330 

MS 2137.3-2353 0.313 1.501 0.778 0.105 

MS 2053.7 0.583 3.146 0.968 0.209 

PKS 0745-191 0.103 0.433 0.818 0.065 

SDSS J0037-0942 0.1955 0.6322 0.6825 0.1026 

SDSS J0216-0813 0.3317 0.5235 0.3632 0.0684 

SDSS J0737+3216 0.3223 0.5812 0.3039 0.0458 

SDSS J0912+0029 0.1642 0.324 0.5325 0.0789 

SDSS J0956+5100 0.2405 0.47 0.414 0.0628 

SDSS J0959+0410 0.126 0.5349 0.5599 0.1152 

SDSS J1250+0523 0.2318 0.795 0.6179 0.0996 

SDSS J1330-0148 0.0808 0.7115 0.7762 0.1184 

SDSS J1402+6321 0.2046 0.4814 0.6575 0.1166 

SDSS J1420+6019 0.0629 0.5352 0.8593 0.1268 

SDSS J1627-0053 0.2076 0.5241 0.5078 0.0779 

SDSS J1630+4520 0.2479 0.7933 0.8114 0.1347 

SDSS J2300+0022 0.2285 0.4635 0.5531 0.0951 

SDSS J2303+1422 0.1553 0.517 0.8651 0.1519 

SDSS J2321-0939 0.0819 0.5324 0.896 0.1312 

Q0047-2808 0.485 3.595 0.8872 0.1606 

CFRS03-1077 0.938 2.941 0.6834 0.1377 

HST 14176 0.81 3.399 0.9757 0.1795 

HST 15433 0.497 2.092 0.929 0.2067 

MG 2016 1.004 3.263 0.5035 0.1234 

SDSS J0029-0055 0.227 0.9313 0.6356 0.1317 

SDSS J0044+0113 0.1196 0.1965 0.3877 0.0573 

SDSS J0109+1500 0.2939 0.5248 0.3803 0.0766 

SDSS J0252+0039 0.2803 0.9818 1.3426 0.2636 

SDSS J0330-0020 0.3507 1.0709 0.8498 0.2109 

SDSS J0405-0455 0.0753 0.8098 1.0851 0.1628 

SDSS J0728+3835 0.2058 0.6877 0.9477 0.1448 

SDSS J0822+2652 0.2414 0.5941 0.6056 0.1004 

SDSS J0841+3824 0.1159 0.6567 0.9671 0.143 

SDSS J0935-0003 0.3475 0.467 0.1926 0.0437 

SDSS J0936+0913 0.1897 0.588 0.6409 0.0953 

SDSS J0946+1006 0.2219 0.6085 0.6927 0.1452 

SDSS J0955+0101 0.1109 0.3159 0.8571 0.159 

SDSS J0959+4416 0.2369 0.5315 0.5599 0.1152 

SDSS J1016+3859 0.1679 0.4394 0.6204 0.0963 

SDSS J1020+1122 0.2822 0.553 0.524 0.0931 

SDSS J1023+4230 0.1912 0.696 0.836 0.1454 

SDSS J1029+0420 0.1045 0.6154 0.7952 0.1231 

SDSS J1032+5322 0.1334 0.329 0.4082 0.0618 

SDSS J1103+5322 0.1582 0.7353 0.9219 0.159 

SDSS J1106+5228 0.0955 0.4069 0.6222 0.0928 

SDSS J1112+0826 0.273 0.6295 0.5052 0.0885 

SDSS J1134+6027 0.1528 0.4742 0.6687 0.1005 

SDSS J1142+1001 0.2218 0.5039 0.6967 0.1735 

SDSS J1143-0144 0.106 0.4019 0.8061 0.1182 

SDSS J1153+4612 0.1797 0.8751 0.7138 0.1305 

SDSS J1204+0358 0.1644 0.6307 0.6381 0.1132 

SDSS J1205+4910 0.215 0.4808 0.5365 0.0803 

SDSS J1213+6708 0.1229 0.6402 0.5783 0.0883 

SDSS J1218+0830 0.135 0.7172 1.0498 0.158 

SDSS J1403+0006 0.1888 0.473 0.6352 0.1332 

SDSS J1416+5136 0.2987 0.8111 0.8259 0.2134 

SDSS J1430+4105 0.285 0.5753 0.509 0.1266 

SDSS J1432+6317 0.123 0.6643 1.1048 0.1662 

SDSS J1436-0000 0.2852 0.8049 0.775 0.1563 

SDSS J1443+0304 0.1338 0.4187 0.6439 0.1 

SDSS J1451-0239 0.1254 0.5203 0.7262 0.1275 

SDSS J1525+3327 0.3583 0.7173 0.6526 0.1611 

SDSS J1531-0105 0.1596 0.7439 0.7628 0.1147 

SDSS J1538+5817 0.1428 0.5312 0.972 0.172 

SDSS J1621+3931 0.2449 0.6021 0.8042 0.1765 

SDSS J1636+4707 0.2282 0.6745 0.7093 0.1276 

SDSS J2238-0754 0.1371 0.7126 1.1248 0.1812 

SDSS J2341+0000 0.186 0.807 1.1669 0.2049 

Q0957+561 0.36 1.41 1.3103 0.1474 

PG1115+080 0.31 1.72 0.7036 0.1604 

MG1549+3047 0.11 1.17 0.5728 0.1194 

Q2237+030 0.04 1.169 0.6685 0.22 

CY2201-3201 0.32 3.9 0.8526 0.305 

B1608+656 0.63 1.39 0.646 0.2154 
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Table 2: Fits to different cosmological models from 80 full D^ s jD s data and 46 restricted 
Dd s /D s data. Fixed value of fl m = 0.27 is assumed except ACDM. 



Cosmological model 


Best-fitting parameters (n = 80) 


Best-fitting parameters (n = 46) 


ACDM 


O - n 20+ 07 




wCDM 


«; = -i-02±g^ 


^ = -1-15^ 


CPL 


wq = 0.60 ± 1.76 


w = -0.24 ± 2.42 




Wa = -7.37 ±8.05 


w a = -6.35 ±9.75 
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W W 

Fig. 5. — The 68.3 and 95.8 % confidence regions for CPL parametrization in the (wo,w a ) 
plane obtained from 80 full D^ s jD s data, 46 restricted D^s/ D s data, and 557 Union2 SNe 
la data. The crosses represent the best-fit points and a star corresponding to ACDM model 
is also added for reference. 



